Absence of diffusion in certain random lattices: Numerical evidence. 



On 
O 
O 
<N 



P. Markos 
Department of Physics, FEI STU, 81299 Bratislava, Slovakia 

Two numerical experiments are performed to demonstrate the physical character of electron local- 
ization in a disordered two-dimensional lattice. In the first experiment we solve the time-dependent 
Schrodinger equation and show that the disorder prevents electron diffusion. Electron becomes 
spatially localized in a specific area of the system. The second experiment analyzes how an elec- 
tron propagates through a disordered sample. In strongly disordered systems, we identify a narrow 
channel through which an electron propagates from one side of the sample to the opposite side. We 
show that this propagation is qualitatively different from the propagation of a classical particle. Our 
numerical analysis confirms that the electron localization is a quantum effect caused by the wave 
character of electron propagation and has no analogy in classical mechanics. 

PACS numbers: 73.23.-b, 71.30.+h, 72.10.-d 
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I. INTRODUCTION 

The electron localization in disordered systems^ is re- 
sponsible for a broad variety of transport phenomena ex- 
perimentally observed in mesoscopic systems: the non- 
Ohmic behavior of electron conductivity, weak localiza- 
tion, universal conductance fluctuations, and strong elec- 
tron localization!^ 

Localization arises in systems with random potential. 
Let us consider the time evolution of a quantum particle 
located at a time t = in a certain small area of the 
sample. For t > 0, the electron wave function scatters 
spatial inhomogeneities (spatial fluctuation of the poten- 
tial) . Multiple reflected components of the wave function 
interfere with each other. As Anderso n^ prov ed, this in- 
terference can abolish the propagation! 4 * 5 * 6 *^ As a result, 
wave function will be non-zero only within a specific area, 
determined by initial electron distribution, and decays 
exponentially as a function of the distance from the cen- 
ter of localization. The probability to find an electron in 
its initial position is non-zero for any time £, even when 
time increases to infinity, t — > oo. 

Similarly to the quantum bound state, the spatial ex- 
tent of the localized wave function is finite. However, the 
physical origin of the localization differs: a bounded par- 
ticle is trapped in the potential well, while localization 
results from interference of various components of the 
wave function scattered by randomly distributed fluctu- 
ations of the potential. 

Localized electrons cannot conduct electric current. 
Consequently, the probability of electron transmission, 
T, through a disordered system decreases exponentially 
as a function of the system length L: T oc exp —L/£. 
The length £ is called localization length. Those materi- 
als that do not conduct electric current due to electron 
localization are called Anderson insulators. 

In spite of significant theoretical effort, our under- 
standing of electron localization is still not complete. 
Rigorous analytical results were obtained only in the 
limit of weak random ness, where perturbation theories 
are applicable! 8 * 9 * 10 * 1 ^ In the localized regime, we do not 



have any small parameter, so no perturbation analysis is 
possible. Also, as we will see below, the transmission of 
electrons is extremely sensitive to the change of sample 
properties. In particular, small local change of random 
potential might cause the change of the transmission am- 
plitude in many orders of magnitude. Clearly, analytical 
description of such systems is difficult. Fortunately, it is 
rather easy to simulate the transport properties numeri- 
cally. In fact, many quantitative data about the electron 
localization was obtained numerically. 7 * 12 * 1 ^ 4 ^ 

In this paper we describe two simple numerical exper- 
iments which demonstrate the key features of quantum 
localization. In Section [IT] we introduce the Anderson 
model that represents the most simple model for study 
of the electron propagation in the two-dimensional sys- 
tem with a random potential. In Sections |III| and [IV| we 



show how randomness influences the ability of an elec- 
tron to propagate at large distances. We solve numeri- 
cally the time-dependent Schrodinger equation and con- 
firm that after a certain time electron diffusion ceases. 
The electron becomes spatially localized in certain part 
of the disordered lattice. This numerical experiment re- 
produces the Anderson's original problemP 

In Section |V| simulate the scattering experiment. We 
consider an electron approaching the disordered system 
from outside, and calculate the amplitude of the trans- 
mission through the sample. Since the transmission de- 
pends on the actual realization of random potential, we 
can, by a small local change of the potential, estimate the 
probability that electron propagates through any given 
sample area. In this way, we investigate the spatial dis- 
tribution of the electron inside disordered sample. For 
weak disorder, we find that electron is homogeneously 
distributed throughout the sample. On the other side, 
in the localized regime we show that electrons propagate 
through narrow spatial channel across the sample. Al- 
though this channel resembles the trajectory of classical 
particles, we argue that the electron still behaves a quan- 
tum particle. We demonstrate the wave character of the 
electron propagation by a simple numerical experiment. 

Both numerical experiments confirm the main feature 
of electron localization: it has its origin in the wave char- 



acter of quantum particle propagation. There is no lo- 
calization phenomena in classical mechanics. 



II. THE MODEL 

Left Figure [T] represents the two-dimensional lattice 
created by regular arrangement of atoms. We consider 
one electron per atom and define its local energy e(r) = 
eo . If the electronic wave functions of neighboring atoms 
overlap, electrons can propagate through the lattice. The 
periodicity of the lattice creates a conductance band, eo — 
4V < E < eo + 4V, 15 where V is given by the overlap of 
electron wave functions located in neighboring sites. 

A disordered two-dimensional lattice is shown in right 
Fig. [I] Now, lattice sites are occupied by different atoms. 
Therefore, both the energy of the electron on a given 
site, e(r), and the hopping term between two neighbor- 
ing atoms, V(r — r'), become position-dependent. In our 
analysis we assume that energies e(r) are randomly dis- 
tributed according to the Box probability distribution, 
P(e) = 1/W if -W/2 < e < W/2, otherwise P(e) = 0. 
We also require these random energies on different sites to 
be statistically independent and assume that the hopping 
amplitude V(r '— r') = V. Although such random lattice 
is rather unrealistic, it imitates all physical features of 
a disordered electron system. The random energies e(r) 
simulate random potential. 

Thus, our random model is characterized by two pa- 
rameters: W represents the strength of the disorder and 
V determines the hopping amplitude. Note that V de- 
fines the energy scale, so we have only one parameter: 
the ratio W/V that we use as a measure of the strength 
of the disorder. 

Let us assume that at a time t = 0a quantum particle 
is located in the position tq. Initial wave function is 



*(r,t = 0) =5(r-r ). 



(1) 



We want to estimate the probability of the electron still 
being in its original position in an infinite time t —> oo. 

The time evolution of the electron wave function is 
determined by the Schrodinger equation, 



ih 



d V{r,t) 

dt 



c(fO*(r,t) + y^;*(P,t), 



(2) 



where \r — r'\ = a is the lattice constant. Equation (|2| 
defines the Anderson model. 

Let us take the zero disorder case, W = first. The 
electron located at time t = in a specific lattice site, 
will diffuse to the neighboring sites. In the limit of an 
infinite time t — > oo, the electron will occupy all sites of 
the lattice. Consequently, the probability to find it in its 
original position equals zero (or, more accurately, it is 
approximately proportional to 1/ (lattice volume)). 

In disordered lattice W ^ 0, electron propagation de- 
pends on the strength of the disorder. Intuitively, one 
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FIG. 1: Left: A regular two-dimensional lattice is a periodic 
arrangement of identical atoms in a rectangular lattice. Right: 
A disordered lattice whose sites are randomly occupied by 
different atoms. The closest distance between two neighboring 
atoms is a. 



expects a very weak disorder not to affect the diffusion 
considerably, but a sufficiently strong disorder should 
stop the diffusion. Then there should be a critical value 
W c : diffusion continues forever when W < W c but ceases 
when W > W c . In the original paper 1 Anderson derived 
the equation for the critical disorder as 



Wc 

V 



= 2eK\n(eK). 



(3) 



According to Eq. (J3|, the critical disorder depends 
only on the lattice connectivity K (the number of near- 
est neighbor sites). Nowadays, we know^iLL^ that the 
dimension d of the lattice is a more important parame- 
ter. In the absence of a magnetic field and of electron 
spin, all states are localized in disordered systems with 
dimension d < d c = 2. Therefore, the critical disorder 
W c = for d = 2 and is non-zero in systems with higher 
dimensionality d > 2. 



III. DIFFUSION 

Now we demonstrate the Anderson's ideas in numeri- 
cal simulation. We examine how the electron diffuses in 
the disordered lattice, defined by Eq. ([2|. The size of 
the system is L x L, where L = 2048a for weakly disor- 
dered samples and L = 1024a for systems with a stronger 
disorder (W/V > 4). 

First, we need to define the initial wave function 
\£(r, t = 0). A more suitable candidate than the 5- 
function (HT) is any eigenfunction of the Hamiltonian de- 
fined on small sub lattice (typically the size of 24a x 24a) 
located in the center of the sampleP^l Usually we chose 
the eigenfunction which corresponds to the eigenenergy 
closest to E = (the middle of the conductance band) . 

To see how the initial wave function develops in time 
t > 0, we solve the Schrodinger equation ([2| numer- 
ically and find the time evolution of the wave func- 
tion \£(r, t). The numerical progr am is based on the 
alternating-direction implicit method^^ used for the so- 
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FIG. 2: (Color online) The quadratic displacement (r 2 (t)) 
(in units a 2 ) as a function of time t. Time is measured in 
h/V. The size of the system is L x L where L = 2048a (L = 
1024a for W/V — 6). Note the logarithmic scale of both axes. 
In weak disorders, we expect the electron to diffuse, so that 
(r 2 (t)) — 2Dt, in accordance with Eq. |5]). Numerically, we 
find that (r 2 ) = 2Dt a with a = 1.004 for disorder W/V = 1 
and a = 0.98 for W/V = 2. The corresponding diffusive 
constants are D = 25.7 and 9.1 (in units a 2 V/K). Only the 
data for time t < 4000H/V were used for W/V = 1, since in 
longer time the electron could reach the edge of the sample. 
The dashed line represents the limit (r 2 ) ma x = L 2 /6, given 
by Eq. (JTl. For stronger disorders, the time evolution of the 
wave function is not diffusive. We find the exponent a ~ 0.82 
(W/V = 4) and a w 0.39 (W/V = 6). 
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FIG. 3: (Color online) The same data as in Fig. 2] but on a 
linear scale. Only the data for small disorder is shown. Note 
that for W/V = 1, (r 2 (£)) is linear only when t < 4O00H/V. 
This is because the electron already reaches the edge of the 
sample. 



lution of elliptic partial differential equations. The algo- 
rithm is described in Appendix A. 

The ability of an electron to diffuse through the sample 



is measured by a quadratic displacement, denned as 

(r 2 (t))=Jdrr 2 \W,t)\ 2 . (4) 

Figures El and p3| show that in weak disorders, W/V = 1 
and 2, \P{t)) is a linear function of time t, 



(r 2 (t))=2Dt. 



(5) 



The parameter D is a diffusive constant which enters the 
Einstein formula for electric conductivity a, 



a = e 2 Dp. 



(6) 



Here e is the electron charge and p is the density of 
states P 

Since we analyze only a lattice of a finite size, we have 
to take into account that the t-dependence of the electron 
wave function might be affected by the finiteness of our 
sample. In this case, we not only observe the diffusion, 
but also the reflection of the electron from the edges. 
Quantitatively, diffusion ([5| is observable only when 



(r 2 W(« (r 2 ) r 



r L/2 r L/2 



= 72 / / (^ + V 2 ) dxdy = 

-L J-L/2 J-L/2 



6 

where (r 2 ) max corresponds to the homogeneously dis- 
tributed wave function, |^(r)| 2 = const = 1/L 2 . 

It might seem that the diffusion of electrons shown 
in Figs. [2] and [3] contradicts the localization theory^ 
that predicts all states to be localized in two-dimensional 
systems. However, this is not the case. The prediction 
of the localization theory concerns the limit of an infinite 
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FIG. 4: (Color online) Quadratic displacement (r 2 (t)) as a 
function of time t/to, to = 1000^/1/ for three systems of the 
size L — 1024a and disorder W/V — 6 (triangles). Although 
(r 2 ) does not increase when time increases, it fluctuates as a 
function of time. The limiting value, R 2 (Eq. |8l depends on 
the actual realization of the random disorder s(r) in the given 



sample, the dashed line shows (r 



L 2 /6 = 174762 a 2 



which is 50 x larger than actual values of (r 2 ). For compar- 
ison, we also show the quadratic displacement for a system 
with stronger disorder, W/V — 8, which is typically 130a 2 . 



system size. Physically, localization occurs only when the 
size of the sample exceeds the localization length, L > £. 
Since £ is very large in weak disorder (£ ~ 10 6 a when 
W = l)}^we observe metallic behavior and diffusion of 
electrons in Fig. [2] Of course, even in the case of W/V = 
1 we would observe localization if much larger systems 
are taken into account P In general, we can observe the 
localization if we either increase the size of the system 
or reduce the localization length. The latter is easier to 
perform, as it requires us only to increase the disorder 
strength W. We will do it in the next Section. 



IV. ABSENCE OF DIFFUSION - 
LOCALIZATION 

The data in Fig. [2] also confirms that the time evolu- 
tion of the wave function is not diffusive when the disor- 
der W increases. Linear increase of (r 2 (t)) is observable 
only for short initial time interval. For any longer time, 
the spatial extent of the electron increases very slowly 
and finally ceases (Fig. [4]). the electron becomes local- 
ized. 

To demonstrate the electron localization more explic- 



itly, we repeat the experiment in Section III with a 
stronger disorder W/V = 6. Similarly to the previous 
experiment, the initial wave function is non-zero in the 
small area 24a x 24a located at the center of the sample. 
For shorter times, we observe that the spatial extent of 
the wave function increases. Then, after a while, (r 2 (t)) 
saturates: 



lim(r 2 (t)) 



R 2 « (r% 



(8) 



Although the spatial distribution of the electron varies 
in time, (r 2 (t)) does not longer increase even if the time 
t increases ten and more times. 

Figures [5] and [6] show the spatial distribution of the 
wave function, |\£(r, t)\. they represent the lattice sites 
with |\£(r)| > 10 -4 . This means that the probability to 
find the electron in any other lattice site is less than 10 -8 . 

Note that there is no potential well in the center of the 
sample where the electron is localized. The only reason 
for the electron being localized in the lattice center is 
the initial wave function, \£(r, t = 0), which was non- 
zero only in the center of the lattice. Applying the initial 
wave function localized in any other area of the sample, 
we would achieve electron localization in that area. This 
is demonstrated in Fig. [7] showing the time development 
of the wave functions of four electrons in the same lattice. 
The initial position of the electrons is centered around 
four points 



x± = L/2 ± L/4, y± = L/2 ± L/4. 



(9) 



We see that in time t > each electron is localized around 
its initial position. This proves that localization is indeed 
the result of interference of wave functions. The electron 
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FIG. 5: (Color online) Spatial distribution of an electron in 
sample with disorder W/V = 6. The size of the lattice is 
1024a x 1024a. Time is given in units of t = lOOOa/V. 
The different colors show sites where |^(r)| > 10 -4 (gray), 
> 5xl0 -4 (brown), 10" 3 (blue), 5xl0 -3 (red), and > 5xl0 -3 
(black). The probability to find an electron on any other site 
is less than 10 -8 . 



is not trapped in any potential well. The localized state 
is not a bound state. 

Figure [7] also shows that the localized states are very 
sensitive to the realization of the random potential. The 
spatial distribution of each electron reflects the local dis- 
tribution of random energies e(f). This is shown quanti- 
tatively in Fig. El where we plot (r 2 (t)) as a function 
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FIG. 6: (Color online) The same as in Fig. [5] only the time is 
t = 500£ and 900t (to = 1000 H/V). 



FIG. 7: (Color online) The time development of four electrons 
located in time t — in four different areas of the same lattice. 
The electrons do not leave the initial areas. The size of the 
sample is L = 1024a. Disorder W/V = 8. Again, time is 
measured in units of to = 1000^/V. 



V. TRANSMISSION THROUGH DISORDERED 

SAMPLE: HOW AN ELECTRON PROPAGATES 

THROUGH DISORDERED SYSTEM? 



of time for three different realizations of the random 
disorder. We see that although all three samples have 
the same macroscopic parameter W/V = 6, the limiting 
value R 2 = lim^^ (r 2 (t)) is not universal but depends 
on the actual distribution of random energies in the given 
sample. Moreover, (r 2 (t)) fluctuates as a function of time 
t. 



Consider now another experiment, frequently used in 
the mesoscopic physics: We take a disordered sample, the 
same as used in previous Sections, to examine what is the 
probability that an electron propagates from one side of 
the sample to the opposite side. Both in experiments and 
in numerical simulations the sample is connected to two 
semi-infinite, disorder- free leads which guide the electron 
propagation towards and out of the sample (Fig. M). An 
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FIG. 8: Schematic description of the scattering experiment for 
measurement of the transmission. The sample is connected 
to two semi-infinite leads represented by regular lattice with 
zero disorder. Inside the sample, the disorder is non-zero. 
If electron comes from the left, it either propagates through 
the sample and contributes to the transmission, or is reflected 
back to the left lead. 



incoming electron either propagates through the sample 
or is reflected back. The probability of transmission, T, 
determines the conductance) 21 * 22 ^ 



9 



-h T ' 



(10) 



Eq. (10) is commonly referred to as the Landauer for- 
mula. It was originally derived for a one-dimensional 
system but can also be used for the analysis of two- and 
more-dimensional samples. Since the width of the leads is 
non-zero the transmission T can be larger than lP^The 
transmission is calculated by the transfer matrix method 
described in RefsF^IMJ 

Contrary to the diffusion problem, discussed in Sects. 
|III|IV[ in the present experiment we do not analyze the 
time development of the electron wave function. Instead, 
we chose the energy E of the electron (E = 0, that is 
the center of the energy band), and calculate the time 
independent current transmission T from the left side of 
the sample to the right side. 

To show how electrons are distributed within the sam- 
ple, we apply Pichard's ideaP^ Let us change the sign of 
a single random energy e(fo) at a site r$: e(fo) — > — e(ro) 
and calculate how this change will influence the total 
transmission T of an electron through the sample. We 
expect that T is sensitive to the change of e(fo) only if the 
electron occupies the site fo, i.e. when |\I/(fo)| is large. 
Contrary, if |^(rb)| is negligible, then the change of e(fo) 
cannot affect the transmission T. 

Thus, by the comparison of the transmissions through 
two systems differing only in the sign of the random en- 
ergy s(ro), we can estimate whether the electron, prop- 
agating through the sample, travels through the site fo 
or not. In repeating this analysis for all lattice sites, we 
can visualize the path of the electron through the sam- 
ple. For numerical reasons, we restricted the sample size 
to 100a x 100a. 

Our results are summarized in Fig. [9j In weak dis- 
orders, W/V = 2, we see that the changing of only one 
random energy has an almost negligible influence on the 
transmission. Typically, T changes only by 1% (or even 
less) when the sign of e(f) changes. Also, all lattice sites 
are more or less equivalent. We conclude that in the 
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FIG. 9: (Color online) Sensitivity of the transmission trough 
the disordered system to the change of the sign of a single 
random energy r - A change of the sign of the random en- 
ergy on orange, red and black sites causes the change of the 
conductance by more than 1%, 10% and 100%, respectively. 
The transmission T is 4.998, 0.52 and 0.00084 for the disor- 
der W/V = 2, 4 and 6 (from top to bottom). The size of the 
system is 100a x 100a, and the electron propagates from the 
left side of the sample to the right side. 



course of the transmission the electron is "everywhere" : 
it propagates through the entire sample as a quantum 
wave. This observation is the key idea of the Dorokhov- 
Mello-Pereyra-Kumar theory of the electron transport in 
weakly disordered systems^ and of the random matrix 
theory of diffusive transport.^ 

The homogeneity of the electron distribution gets lost 
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FIG. 10: (Color online) The same as in Fig. |9]but disorder 
W/V = 10. The transmission T = 9 x 10~ 15 . 
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FIG. 12: (Color online) The electron path through two 
strongly disordered samples: both samples have the same real- 
ization of random energies. They differ only in the amplitude 
of fluctuations. Shown are the lattice sites where the change 
of the sign of random energy causes the change of the loga- 
rithm of the transmission by 1% (orange) and 10% (red). We 
see that the electron prefers completely different trajectories 
through these samples. 
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FIG. 11: (Color online) The same as in Fig. [9] but for disor- 
der W/V — 20. The transmission is really small, InT = —96. 
The change of the sign of the random energy on gray, brown, 
orange and red sites causes a change of the logarithm of the 
transmission by more than 0.01%, 0.1%, 1% and 10% respec- 
tively. Although it seems that the path through the sample 
is determined by a valley in the potential landscape, this is 
not the case. The inset shows sites of the sample where the 
random energy |e| < 1. 



as the sites on the main trajectory (Fig. 10). This in- 



dicates that the electron still feels the entire sample and 
its propagation is highly sensitive to any change of the 
realization of the random potential. 

the resulting trajectory cannot be identified with any 
valley or equipotential line in the random potential land- 
scape. To demonstrate this, we show in Fig. [TTJthe tra- 
jectory of an electron through an extremely strongly dis- 
ordered system (W/V = 20 - in this case, we consider the 
change of the logarithm of the conductance) . Although 
the trajectory of an electron seems to be well defined, 
there is no continuous potential valley which might sup- 
port the propagation. The inset of Thus, the choice of 
the transmission path is the result of quantum interfer- 
ence: an electron arising from the left inspects the entire 
sample and finds the most convenient spatial "channel" 
for its propagation. We cannot speak about a trajectory 
in the sense of classical particles. 

To support our last claim, let us consider two samples, 
with the same realization of random energies e(r), but 
different amplitudes of random energies: W/V = 10 for 
the sample I and W/V = 20 for the sample II: 



when the disorder increases! 26 * 27 * The change of the ran- 
dom energy sign on some sites influences the transmission 
more than the same change on other sites. Some areas of 
the sample seem not to be visited at all. We can see the 
formation of the electron "path" through the sampled 
This path is clearly visible for very strong disorder shown 
in Figs. [10) [IT] and [12] 

However, we want to stress that even in case of strong 
disorders we cannot speak about the path in its classi- 
cal sense. Even if the electron path is well visible, there 
are still other sites, often located on the opposite side of 
the sample, that influence the transmission as strongly 



e(r) n = 2e(r) 1 



(11) 



for each site r. With the help of the above-mentioned 
method, we find the trajectories of electrons through 
these two samples. For the propagation of classical par- 
ticle both trajectories (for the sample I and sample II) 
should coincide. However, an electron is not a classical 
particle. As shown in Fig. [12] the ways how an elec- 
tron propagates through the two samples, I and II, are 
completely different. An increase of fluctuations of the 
random potential causes the electron to choose a com- 
pletely different route. 



VI. CONCLUSION 

We discussed two features of localization of quantum 
particle in a disordered sample. Firstly, we demonstrated 
numerically that the diffusion of the quantum particle 
through randomly fluctuating potential ceases after cer- 
tain time. The particle becomes spatially localized. The 
physical origin of localization is different from the bound- 
ing of a particle in a potential well. Localization is caused 
by a multiple scattering of the wave function on randomly 
distributed impurities (fluctuations of the random poten- 
tial) . It is not due to the trapping of the particle in the 
potential well. 

In the second part of the paper, we examined the prop- 
agation of a quantum particle through a disordered sam- 
ple and discussed how this propagation depends on the 
disorder. Again, we confirmed indirectly the wave char- 
acter of the propagation. This drove us to the conclusion 
that the electron localization is a purely quantum effect 
without any analogy in classical mechanics. 

In both numerical experiments, the key condition for 
the localization to happen is the quantum coherence of 
the wave function. This is generally not fulfilled in ex- 
periment, where the incoherent scattering - for instance 
the scattering of electrons with phonons - plays the cru- 
cial role. As any incoherent scattering destroys quantum 
coherence, the observing of electron localization experi- 
mentally requires the mean free path of incoherent scat- 
tering to be larger or at least comparable to size of the 
sample. This happens at a very low temperature. Of 
course, localization does affect the transport of electrons 
also at higher temperatures. These effects are, however, 
above the scope of present discussion. 

With localization being a wave phenomenon, we can 
expect the similarity of quantum propagation with clas- 
sical wave phenomena. 29 That enables us to observe lo- 
calization in many other instances. In particular, we can 
expect that classical waves - electromagnetic or acous- 
tic - will also be localized in a disordered mediumP^ 
The localization of microwave electromagnetic waves 
was experimentally observed 31 . Another very interest- 
ing experiment 32 proves the weak localization of seismic 
waves. 



APPENDIX A: NUMERICAL SOLUTION OF 
THE SCHRODINGER EQUATION 

To integrate the Schrodinger equation Q numerically, 
we first have to discretize the time derivative: 



9*(r,t) 



[^(r,t + St)-^(r,t)}, (Al) 



dt St 

where St is the time iteration step. In order to define 
the iteration procedure for the calculation of the wave 
function in \£(r, t + S) in terms of \£(r, £), we also have 
to determine the time in which the wave function on the 
right-hand side of Eq. pi is calculated. The explicit 
method takes the entire r.h.s. of Eq. Q in time t. 



The resulting iteration scheme is simple, the time step 
St, however, has to be very small to avoid numerical 
instabilities.^ Other, more sophisticated, explicit meth- 
ods based on the Suzuki- Trotter formula are described 
in Ref. 18 In this paper, w e use the alternating-direction 
implicit iteration schema .H212QJ This method presses the 
numerical integration of Schrodinger equation ^ in two 
steps. Firstly, we write the Schrodinger equation pi) in 
its discrete form (r = (x,y)) 

^{x, y, t + St) = #(*, y) + —St [ye(x, y)V(x, y, t) 

+^(x + a,y,t + St) + Stf(x -a,y,t + St) 



H-^0, y + a, t) + #(a, y-a,t) 

(A2) 

Note that the wave function along the x direction is in 
time t + St, while the wave function in y direction is in 
time t. The second step is to put t — > t + St and re- write 
the iteration equation in the form 

tf (a, y, t + St) = tf (a, y) + —St [ye(x, y)V(x, y, t) 

+Sfr(x + a, y, t) + V(x - a, y, t) 

+^(x, y + a,t + 6t) + ^(x, y - a, t + St) 

(A3) 
Now the wave function along the y direction is given in 
time t + St and the wave function along x direction is in 
time t. 

The advantage, that this algorithm presents is its nu- 
merical stability, even for rather large values of St. The 
price for the numerical stability is that we have to solve 
N (N = L/a) systems of N linear equations in each it- 
eration step. Fortunately, iteration schema ( A2|A3 ) only 
requires solution of the three-diagonal system of linear 
equations, what is easy to calculate. 

Other modifications of the above algorithm are pos- 
sible since we are free to sub stitute the wave function 
^f(x,y,t) on the r.h.s of Eqs. ( |A2|A3[ ) by ^(x,y,t + St), 
or, eventually, by \^{x,y,t) + ^(x,y,t + St)]/ 2. 

To measure the accuracy of the numerical solution, we 
also calculate in each time t the norm of the wave func- 
tion, 



M-- 



/ df|*(r, 



X 4\|2 



t)\\ 



(A4) 



We obtained that for St = O.lh/V the norm JV is always 
very close to 1: 



\J\f-l\ <2x 1(T 3 . 



(A5) 
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